Ergodicity from Nonergodicity in Quantum Correlations of Low-dimensional Spin 

Systems 
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Correlations between the parts of a many-body system, and its time dynamics, lie at the heart of 
sciences, and they can be classical as well as quantum. Quantum correlations are traditionally viewed 
as constituted out of classical correlations and magnetizations. While that of course remains so, we 
show that quantum correlations can have statistical mechanical properties like ergodicity, which is 
not inherited from the corresponding classical correlations and magnetizations, for the transverse 
anisotropic quantum XY model in one-, two-, and quasi two-dimension, for suitably chosen trans- 
verse fields and temperatures. The results have the potential for applications in decoherence effects 
in realizable quantum computers. 



I. INTRODUCTION 

Quantum entanglement [1], a striking property of 
quantum theory, plays a crucial role in many quantum in- 
formation processes including achieving secure quantum 
cryptography [2||, faster quantum computers, and bet- 
ter quantum communication [3j . Quantum entanglement 
in many-body systems has created considerable interest 
recently because of its potential in understanding quan- 
tum many body phenomena, vifhich are not amenable to 
classical perceptions, such as quantum phase transitions 
[1], decoherence Q, etc. Also, quantum entanglement 
forms a central element in bringing together many-body 
physics and quantum information science, and hence in 
the development of new exciting physics [B-Q- Prop- 
erties of entanglement in many-body systems like cold 
atomic gases in optical lattices [9] , trapped gaseous Bose- 
Einstein condensates ilfl,] , and spin models [1, 0, ITTl - flsj 
have been dealt with increasing interest in recent times. 

The importance of studying the time-evolved state can- 
not be over-emphasized in both many-body physics and 
quantum information. In particular, important processes 
in quantum information, like the one-way quantum com- 
puter crucially depends on the time-evolved state and the 
possibility to control it in experimental set-ups [13] ■ On 
the other hand, the changes happening in several physical 
parameters, including statistical mechanical properties of 
such parameters, during evolution, is an important as- 
pect in many-body physics, to understand, for example, 
the phenomena of decoherence in the system How- 
ever, despite its vital importance, the complexity of the 
time-dynamics in physically interesting models, has led 
to a limited amount of study of such phenomena 0, 0] ■ 

The validity of a statistical mechanical description of a 
physical quantity depends on the behavior of that quan- 
tity in the time-evolution of the system. Ergodicity is a 
necessary condition for the validity of a statistical me- 
chanical description of a physical quantity. A physical 
property is said to be ergodic if the time average of the 
quantity matches its ensemble average. 

A two-party quantum state (pure or mixed) has non- 
zero classical correlations whenever it has non-zero quan- 
tum correlation with respect to any "good" entanglement 



measure [15| (cf. [l6|). However we show that the prop- 
erties of quantum correlations may not be traced back to 
those in classical correlations. In particular, in the case 
of bipartite quantum states, which in our case are near- 
est neighbor states of quantum spin systems, we show 
that statistical mechanical properties like ergodicity can 
appear in quantum correlation even when it is absent for 
all classical correlations and magnetizations. 

To deal with questions about ergodicity, we consider 
the anisotropic transverse quantum XY model in one- 
dimension, ladder, and two-dimension. Spectacular ad- 
vances in cold gas experimental techniques have made 
these models experimentally realizable in such systems 
[1, Q . The two-component Bose-Bose and Fermi- Fermi 
mixture, in the strong coupling limit with suitable tuning 
of scattering length and additional tunneling in the sys- 
tem can be described by the quantum XY Hamiltonian. 
The dynamics of the system can be simulated by control- 
ling the system parameters and the applied transverse 
field. Therefore the phenomenon discussed in the paper 
can potentially be verified in the laboratory. Moreover, 
many solid state compounds can be described well by 
this model. In particular, the LiHo2:Yi_a;F4 compound 
is known to be described by the three-dimensional (quan- 
tum) transverse Ising model, which is a specific case of 
the XY model 

The question of ergodicity in the quantum XY chain 
was already considered by Mazur in Ref. 18] for magneti- 
zation, and important further developments were carried 
outinRcfs. [la, [13 for magnetization as well as classical 
correlations, and all these physical quantities were shown 
to be nonergodic. Here we show that the correspond- 
ing quantum correlation aka entanglement can however 
be ergodic. Moreover, we carry over the considerations 
to other low-dimensional spin systems on ladders and in 
two-dimensions, and in all cases we show that quantum 
correlation can be ergodic without the classical correla- 
tions and magnetizations being so. 

Quantum correlations are known to provide us with ap- 
plications in quantum information processing tasks, like 
quantum communication, quantum cryptography, and 
quantum computation, that cannot be matched by sys- 
tems with classical correlations only [l|-[ll- Also, some 
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form of quantum correlation can potentially be a "uni- 
versal" detector of quantum phase transitions 0] , and 
numerical algorithms inspired by quantum information 
concepts can potentially provide better understanding of 
many-body systems [21.] . And yet, statistical mechanical 
properties like quantum phase transitions can be identi- 
fied by using both classical |3| and quantum correlations 
[g!, 7]. However, the fact that even statistical mechanical 
properties like ergodicity can emerge in quantum corre- 
lation of a quantum state without it being present in 
the classical correlations (and magnetization) had so far 
eluded us. The findings can potentially be important 
in several areas in quantum many-body physics and in 
quantum information. In particular, it may turn out to 
be crucial for studies on the decoherence times in quan- 
tum computers [22 1. 



II. THE XY MODEL IN THE TRANSVERSE 
FIELD 

A. Description of the model 

We consider an interacting spin Hamiltonian with an 
external magnetic field given by 



H — Hi, 



h{t)H, 



mag J 



where h{t) is a time-dependent function. To obtain non- 
trivial effects on the evolution due to the external mag- 
netic field, one must choose the magnetic field and other 
parameters in such a way that the interaction part and 
the field part of the Hamiltonian do not commute, i.e., 
[iJint, -ffmag] 7^ 0. A simple way to obtain this is to choose 
the field part as 



and the interaction term as 



The paper is organized as follows. In Sec. |TI1 we give 
a brief description of the considered model and the time- 
dependent transverse field. In one-dimension (ID), this 
model can be exactly diagonalized jl^, [2^ [1^ , and we 
perform the corresponding Jordan- Wigner transforma- 
tion to arrive at the expressions for magnetization, clas- 
sical correlations, and entanglement for thermal (canon- 
ical) equilibrium states as well as for time-evolved states 
for an infinite chain. In this paper, we will use logarith- 
mic negativity as the measure of entanglement, and we 
define this in Sec. IIIII In Sec. IIVI we consider the ques- 
tion of ergodicity in the infinite chain, and we show that 
there are situations where magnetization and all classi- 
cal correlations are nonergodic, while entanglement re- 
mains ergodic. Additionally, we show that in a large 
number of the remaining cases, entanglement is ergodic 
while magnetization and almost all classical correlations 
are nonergodic. Entanglement, therefore, follows the mi- 
nority verdict rather than the majority one: Although 
most of the classical correlations and magnetization are 
nonergodic, entanglement is actually ergodic. In Sec. |Vl 
we show that the statistical behavior of entanglement, 
classical correlations, and magnetization for finite XY 
chains is similar to that of the infinite ID XY model, 
even for relatively moderate system size. Specifically, we 
consider finite chains of size up to iV = 12, with periodic 
boundary conditions. The qualitative similarity of the 
results obtained for finite and infinite chains led us to 
consider the statistical properties of finite spin systems 
in two-dimensional (2D) and quasi-2D systems with pe- 
riodic boundary conditions, where an exact diagonaliza- 
tion procedure is not possible in the infinite limit. We 
find that the ergodicity of entanglement, despite magne- 
tization and all classical correlations being nonergodic, 
persists in higher dimensions (Sec. I VII and Sec. IVIip . We 
summarize our results in the final section (Sec. IVIII[) . 



jY.[{i + i)sts^^+{i~i)sjsy] 



with 7 7^ 0, where J is the coupling constant, 7 is the 
anisotropy measure, Sf , Sf and are one- half of the 
Pauli spin matrices at the i*^ site, and (ij) denotes inter- 
actions between all the nearest-neighbor pairs in a lattice. 
We assume periodic boundary conditions. Therefore the 
time-dependent total Hamiltonian that we study in this 
paper, reads 

H = jJ2[il+ l)S'tS] + (1 - 1)8^3^] - h{t) St, 

(1) 

and is called the anisotropic (quantum) XY model in a 
transverse field. 

The transverse field applied is in the form of an initial 
disturbance that subsequently goes to zero for t ^ 0: 



hit) = { I 



t < 
t > 



(2) 



where a ^ 0. Hence we study the system that has evolved 
with time after an initial disturbance due to the field. 

The state that we consider evolves according to the 
Hamiltonian H, given in Eq. ([IJ. But it also depends 
on the initial state, from which it starts evolving. We 
assume that the initial state of the evolution is the (ther- 
mal) equilibrium state at the initial time, and at absolute 
temperature T. Let us denote the canonical equilibrium 
state at time t as Peq(t), so that 



exp[-/3i/(t)] 



Here Z is the partition function, given by 

Z = ti[exp[- /3H{t)]], 
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and f3 = where k is the Boltzmann constant. To 
consider questions about ergodicity, we will be interested 
in the behavior of the evolved state. The evolution is 
governed by the Hamiltonian H{t), from a given initial 
state. Since we will compare the long-time average of the 
properties of the evolved state, with those of the equilib- 
rium state, it is natural to suppose that the initial state 
is the equilibrium state at t = 0, so that there is a possi- 
bility that the physical quantity after a sufHciently long 
time, will equilibriate. We denote the evolved state by 
where the prefix corresponds to the temperature 
of the initial equilibrium state ^"^(0) — p"(0). 

B. Single- and two-site reduced density matrices 

The main intention in this paper is to study the statis- 
tical behavior of single-site magnetization, nearest neigh- 
bor entanglement and correlations in various spin models. 
Therefore we now find out the single-site and two-site re- 
duced density matrices of the equilibrium and evolved 
states. 

A general single-site density matrix is given by 

p'^ll + 2A/ • S, 

where / denotes the 2x2 unit matrix and magnetization 
is given as M = tr[p^tT], where a = {a^,a'^,a^) are the 
Pauli matrices. The single-site density matrix can be 
obtained by tracing out N — 1 sites from the TV-site spin 
states. Let us now consider the single-site density matrix 
for the equilibrium state p^q{t). Due to symmetry, all the 
single-site density matrices of the equilibrium states are 
equal. We will denote them by piq{t) (hiding the prefix 
/?). Now piq{t) — piq{t), where the complex conjugation 
is taken in the computational basis, which (for each site) 
is the eigenbasis of the Pauli matrix a^. Therefore 

Myq{t)^tT[sypiq{t)] = o. 

Moreover the Hamiltonian H{t) has the global phase flip 
symmetry {[H,IliS^] =0), from which it follows that 

M:^{t)^tT[s'^pi^{t)] = o. 

Consequently, the single-site density matrix of the equi- 
librium state reduces to 

piq{t) = ^I + 2M^q{t)S\ 

The evolved state does not necessarily have the property 
of being equal to its complex conjugation, and consider- 
ation of the global phase flip symmetry is complicated 
by the fact that the Hamiltonian is explicitly depen- 
dent on time. However, using the Wick's theorem, as 
in [l^, [23, HI] , the single-site density of the evolved state 
turns out to be of the form 

p'it) = ^I + 2M'{t)S'. 



So, the single-site (transverse) magnetization of the equi- 
librium state is 

M4(t)-tr[5VeV0], 

while that for the evolved state is 

M'^it) ^tT[S'p^it)]. 

The general two-site density matrix for the eqilibrium 
and the evolved state can therefore be written as 

where T*-' — tr[p-'^((T* igj a^)] represents two-site corre- 
lation functions. Using Wick's theorem, once again we 
can show that the yz and xz correlations will vanish for 
the evolved state. For the equilibrium state, only the 
diagonal correlations, T^^, T^^, and T^^, remain. 

C. Two-site reduced density matrix for the ID 
infinite spin model 

The Hamiltonian for the ID infinite XY spin model 
with nearest neighbor interaction is given by 

^ = I E [i^+^y^<'^+l + (1 - 7)'^r^r+l -Ht)at] . (3) 

i 

This model is exactly solvable by successive Jordan- 
Wigner, Fourier, and Bogoliubov transformations [23j . 
Below we obtain the two-site density matrices for the 
evolved and equilibrium states of the above Hamiltonian 
H. 



1. Evolved state 

Suppose that the quantum system described by the 
Hamiltonian H starts off from the initial state which is 
a (canonical) equilibrium state at temperature T. We 
are interested in the nearest-neighbor (two-site) density 
matrix of the evolved state at time t — )> 00, that started 
off from the equilibrium state. 

The correlations and transverse magnetization in this 
case are given as follows [IB, i20|, [l^l • 

= G(i), 

T"" = [M^'Y ~ G{l)G{~l) 

where G{R) (for R = ±1) is given by 

^^Dx 1 r j^tanh(^A(5)/2 , . . , 2 - n 

G{R) = - / # Ar\A2(n\ ^ sm 0i? sm cj, - cos^ cj,) 
ttJo A(a)A^(0) 

X (7^^ sin^ -I- (cos — a) cos 0) . (5) 
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The transverse magnetization is given by 



A(a)A2(0) 



X COS (/<[(cos — a) cos + 7^ sin^ 



Here 



and 



21 2 



A(x) = {7^ sin^ (/) + [a; - cos(/)]^} 



a=-, /3 = /3J. 
Note that a and /? are dimensionless variables. 



(6) 
(7) 
(8) 



2. Thermal state 

The correlations and transverse magnetization, for the 
thermal equilibrium state at a temperature T, are given 
as follows [Tgl [20l . [23} . The only correlations that remain 
are the diagonal ones. 

r--(a,/3) = G(-l,a,/3), 

Tyy{a,P)^G{l,a,P), (9) 
r-(a, /3) = [M-(a, - G(l, a, /3)G(-1, a, 

where G{R, a, (3) (for i? = ±1) are given by 
1 r , ,tanh(/3A(a)/2) 



G{R,a,l3) 



TT Jo A(5) 
X (7 sin(0i?) sin (/> — cos ^(cos </) — d) ) 



(10) 



and 



M'{a,(3)) 



1 r tanh(/3A(S)/2)(coS(?!)-a) 



(11) 



III. MEASURE OF ENTANGLEMENT: 
LOGARITHMIC NEGATIVITY 

The bipartite entanglement measure that we will con- 
sider throughout the paper is the logarithmic negativity 
(LN) [24l|. It should be stressed, however, that the re- 
sults do not depend on the choice of the entanglement 
measure. To define logarithmic negativity, let us first in- 
troduce negativity. The negativity N{pab) of a bipartite 
state pAB is defined as the absolute value of the sum of 
the negative eigenvalues of , where p^ denotes the 
partial transpose of pab with respect to the A-pait [25| . 
The logarithmic negativity is defined as 



En (pab) = log2[2N (pab) + 1]- 



(12) 



In our case, the bipartite states are states of two qubits, 
so that has at most one negative eigenvalue j26j . 
Moreover, for two-qubit states, a positive LN implies that 
the state is entangled and distillable l25l [27| . while En = 
implies that the state is separable |25| . 



IV. INFINITE QUANTUM XY SPIN CHAIN IN 
A TRANSVERSE FIELD 



In this section, we deal with the question of ergod- 
icity of entanglement and other physical quantities for 
the (one-dimensional) infinite XY spin chain in the 
time-dependent transverse magnetic field, for which the 
Hamiltonian is given by Eq. ([T|), to be considered for 
ID here. The Hamiltonian depends on several physical 
parameters: the coupling J, the anisotropy 7, and the 
initial transverse magnetic field a. To check for ergodic- 
ity of a given physical quantity Q, we begin by looking at 
its behavior in the time-evolution of the system, given by 
a set of system parameters (J, 7, a), and where the initial 
state of the evolution is the canonical equilibrium state of 
the system at the initial time and at a given temperature 
T. We are interested in the time-averaged value of Q 
at large times. Let us call it Q°°{T,a), keeping in mind 
that it also depends on the system parameters J and 
7. The physical quantity Q will be called nonergodic for 
the system under consideration, if its large-time value for 
the canonical states, Q'^^'^{T' ,h{t — 00)), corresponding 
to any value of temperature, T', in some relevant range 
of temperature around T, does not match Q°°{T,a). In 
other words, 

g°°(r, a) + g'='^"(T', h{i = oo)) vt', (13) 

where T' is chosen in a physically relevant range, will 
imply nonergodicity of Q. Otherwise, it will be termed 
ergodic. We term a physical quantity as strongly noner- 
godic if the relation holds for all temperatures T'. 

We begin by looking at the (transverse) magnetization 
in this model. The results and discussions, along with 
the figures presented, are for 7 = 1/2. However, all the 
results are independent of the value of 7. Choosing the 
temperature of the initial canonical state at an exemplary 
value, given by /3 = 20, we find the long-time average of 
magnetization in the time-evolved state of the system for 
different values of the initial magnetic field a. As exem- 
plary cases, for a/J — 0.2, 0.6, and 2, the long-time aver- 
ages of the transverse magnetizations in the time-evolved 
states, where the dynamics is assumed to start off from 
the canonical equilibrium state with /? — 20, are 0.005, 
0.079, and 0.643. Actually, these long-time averages of 
magnetization are all nonzero, as long as the initial mag- 
netic field is nonzero. The canonical magnetization (i.e., 
the magnetization in the canonical equilibrium state) is 
however zero for all nonzero times, and at all tempera- 
tures, impl ying th at the magnetization is nonergodic in 
this model |18l - [2(j |. The difference between the magneti- 



5 



zations of the evolved and the equihbrium states directly 
depends on the applied field strength. 

The situation for the nearest neighbor correlations, 
rpxx ^ rpyy ^ g^j^j j^zz ^ ^ more involved. For low values 

of the initial magnetic field, the correlations remain er- 
godic. As we crank up the initial transverse field, the yy 
and zz fields behave similarly, and become strongly non- 
ergodic, just like the transverse magnetization: There is 
no temperature for which the corresponding equilibrium 
values match the long-time averages. The situation for 
the yy correlation is shown in Fig. 1 (middle). The xx 
correlation is ergodic for low magnetic fields. However, 
for high magnetic fields, T^^ attains ergodicity only if 
we allow temperatures much higher as compared to the 
initial temperature. See Fig. 1 (top). For example, for 
a/ J — 2, the intersection (of the equilibrium curve and 
the time average horizontal line) is at /? < 0.5, which is 
two orders of magnitude below the inital /3 = 20. This is 
not physically meaningful. Therefore, if we suppose that 
we allow the temperature of the equilibrium state to be 
within about one order of magnitude around the initial 
temperature, the equilibrium and time-evolved curves for 
T^^ will not intersect, and we obtain nonergodicity of the 
XX correlation in this model. 

Surprisingly, the situation is the opposite for quantum 
correlations [2^ [2^. The long-time average of the nearest 
neighbor entanglement in the time-evolved state (with 
the initial state of evolution being the canonical state 
with the same temperature as in the case of correlations 
and magnetization) is decreasing with increasing initial 
magnetic field. And for example, all the three horizontal 
lines, in Fig. 1 (bottom), for time-averaged entangle- 
ment corresponding to a/ J — 0.2, 0.6, and 2, intersect 
the curve corresponding to the canonical entanglement 
at different temperatures (within one order of magnitude 
of the initial temperature), implying that entanglement 
is actually ergodic in these models. 

Ergodicity can therefore appear in a physical quantity 
(entanglement, here), without it being present in the con- 
stituent physical quantities (the correlations and magne- 
tization, here). 



V. FINITE XY SPIN CHAIN IN A 
TRANSVERSE FIELD 

Having found the status of the ergodicity of magne- 
tization, and classical and quantum correlations for the 
transverse XY chain, we would like to look at its charac- 
ter for higher dimensional systems. In these cases, how- 
ever, there are no analytical diagonalizations possible, 
and we have to resort to numerical or approximate meth- 
ods. We will be using exact diagonalizations for finite pe- 
riodic higher dimensional systems, to find the status of 
ergodicity of magnetization, and classical and quantum 
correlations in such systems. But before that, in this sec- 
tion, we consider finite periodic chains, and try to check 
whether the inferences on ergodicity remains the same as 
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FIG. 1. (Color online) Correlations (dimensionless) and en- 
tanglement (in ebits) for the infinite ID XY model. The hor- 
izontal axis represents the dimensionless parameter /3. All 
the classical correlations are ergodic at low fields. The cor- 
relations T^^ (middle) and T^^ become nonergodic at high 
initial fields. The xx correlation (top) is ergodic at low fields, 
while at high fields, it is nonergodic if we assume that the 
equilibrium temperature can only be an order of magnitude 
different from the initial temperature. The opposite is true 
for entanglement. The curves with black circles represents the 
equilibrium values as functions of /3J. For easy comparison, 
we have plotted the horizontal lines to show the values of the 
corresponding physical parameter for different values of the 
initial magnetic field; precisely, for a/ J = 0.2 (red squares), 
a/ J = 0.6 (blue diamonds), and a/ J = 2 (pink triangles). 
a/ J is dimensionless. 

those from infinite chains. 

The Hamiltonian for a (one-dimensional) spin chain of 
12 spins, with periodic boundary conditions, is given by 

12 12 

H^JJ2[{1 + i)SfSf+, + (1 - j)SfSf^,] ~h{t) J2 St. 

i=l i=l 

(14) 

Here, we have assumed periodic boundary conditions, 
whereby 5*12+1 = Si. 

As in the case of the infinite chain, we start the evolu- 
tion with the canonical equilibrium state at a tempera- 
ture given by /3 = 20, and compare the long-time average 
of a certain physical quantity with the values of the same 
quanitity in the canonical equilibrium state at different 
temperatures. 

As before, the picture is simple for magnetization. It 
is always zero for the equilibrium state at any nonzero 



time, while it has a nonzero time-averaged value in 
the time-evolved state for any nonzero initial transverse 
field. Hence, magnetization is strongly nonergodic for all 
nonzero a. The behavior of all the classical correlations in 
the finite chain is also similar to that in the infinite chain. 
Of the classical correlations, let us, for definiteness, con- 
sider the yy correlation, T^^. As seen in Fig. 2 (left), it 
remains ergodic for low values of the applied field, while 
becoming strongly nonergodic for higher values. How- 
ever, Fig. 2 (right) shows that with increasing applied 
magnetic field, the entanglement becomes ergodic. 

It is clear therefore that already at the level of TV = 12 
spins of the one-dimensional XY model, the statistical 
mechanical behavior of nearest neighbor entanglement 
and classical correlations, and single-site magnetization, 
mimics those of the infinite chain. This convinces us to 
study such statistical behavior in the case of similar sys- 
tems in higher dimensions, where choosing a finite num- 
ber of particles is a necessity. This will be taken up in 
the two succeeding sections. 




FIG. 2. (Color online) Correlation (dimensionless) and en- 
tanglement (in ebits) for a chain of 12 spins, with periodic 
boundary conditions. The depictions are as before. Clearly, 
the correlation T^^ is ergodic for low magnetic fields, while 
becoming nonergodic for higher values of the field. The op- 
posite is vafid for entanglement. 



VI. QUASI-2D: XY MODEL WITH A 
TRANSVERSE FIELD ON A LADDER 

A ladder of spins forms an interesting intermediate 
dimension between one and two dimensions, containing 
several intriguing features of both the worlds. Moreover, 
there are materials that have been discovered which can 
be described by a ladder Hamiltonian. Currently avail- 
able techniques in cold gas systems also offer the possi- 
bility of realizing such models [s^l . 

Eight spins arranged in the form of a ladder is shown 
in Fig. [3l where we have assumed periodic boundary 
conditions along the rails, so that and 5*1, and Ss and 
are connected by quantum XY interactions just like 
all other nearest neighbor spins. Let us now study the 
statistical mechanical behavior of nearest neighbor corre- 
lations and entanglement, and single-site magnetization, 
for such a spin model. In Fig. SI we plot the nearest 
neighbor entanglement and the classical correlation T^^, 
for this system, as a function of temperature for various 




FIG. 3. (Color online) A ladder of 8 spins with quantum 
XY rails and rungs. Periodic boundary condition is assumed 
along the rails. 
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FIG. 4. (Color onUne) Correlation (dimensionless) and entan- 
glement (in ebits) on a ladder with periodic boundary condi- 
tions. Green stars represent a/J — 1.2. The other depictions 
are the same as before. At high magnetic fields, correlations 
become nonergodic, while entanglement is ergodic. 



applied transverse field strengths, and where again we 
assume that the evolution starts off from a temperature 
given by /3 = 20. 

Magnetization: The transverse magnetization is clearly 
nonergodic, as the equilibrium value is zero for all 
temperatures, while the long-time averaged values are 
nonzero for any nonzero applied initial field. 

Correlations: For the ladder of 8 spins, let us begin by 
considering the behavior of the T^^. See Fig. |4] (left). 
The long-time canonical equilibrium T^^, is vanishing at 
zero (3, which then decreases to a value of —0.109, and 
converges to —0.088 for low temperatures. Now the long- 
time average value of the evolved state T^^, for an applied 
transverse field a/ J = 0.2 is —0.094, which is represented 
by an horizontal line (of red squares) in Fig. |4] (left). It 
clearly indicates that there exists a temperature for which 
T^y of the canonical equilibrium state coincides with the 
long-time average of T''*' of the evolved state, predict- 
ing ergodicity of T*'*'. However, for higher values of the 
applied field strength, the situation is the opposite, and 
for example, for a/J — 0.6 and a/ J — 1.2, the respective 
horizontal lines do not intersect the equilibrium curve, in- 
dicating strong nonergodicity of the T^^ correlation for 
these models. For an even higher field oi a/J — 2, the 
intersection is at a temperature that is more than an or- 
der of magnitude different from the initial temperature, 
and we predict nonergodicity again. The correlation T^^ 
is ergodic at low fields, and strongly nonergodic for high 
fields. The correlation T^^ is again ergodic at low fields, 
but at high fields, it is nonergodic. 

Entanglement: The behavior of entanglement is ex- 
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actly the opposite to that of classical correlations. The 
long-time average of entanglement decreases with in- 
creasing magnetic field. As seen in Fig. |4] (right), for the 
exemplary values, a/J = 0.2,0.6,1.2 and 2, of the initial 
magnetic field, the behavior of entanglement is clearly 
ergodic. 

Thus, it is clear that for the spin ladder (quasi-2D) 
model with quantum XY rails and rungs, the statisti- 
cal behavior of nearest neighbor entanglement and clas- 
sical correlations, and single-site magnetization, is qual- 
itatively the same as that of the spin chain: Quantum 
correlations can be ergodic without classical correlations 
and magnetization being so. 



En 




FIG. 6. (Color online) Correlation (dimensionless) and en- 
tanglement (in ebits) on a 2D square lattice. The depictions 
are the same as before. And as before, high magnetic fields 
render the correlations as nonergodic, while entanglement is 
ergodic there. 



VII. TRANSVERSE XY MODEL IN 2D 

The quantum XY model in a transverse field in two di- 
mensions is an important system in several applications. 
In particular, in the special case of unit anisotropy, the 
dynamics of this model is utilized in the measurement- 
based model of quantum computation • 




FIG. 5. (Color online) A two-dimensional array of 12 spins 
with XY interactions. Periodic boundary conditions is as- 
sumed, so that the system forms a torus. 

We consider a two-dimensional array of 12 spin-1/2 
particles arranged on a square lattice, and interacting 
via nearest neighbor interactions. A schematic represen- 
tation is given in Fig. [H where periodic boundary con- 
dition is assumed, so that the system forms a torus. We 
now study the ergodicity properties of this torus for near- 
est neighbor entanglement and classical correlations, and 
single-site magnetization. In Fig. [SJ we plot the nearest 
neighbor entanglement and the T^^ correlation for the 
torus of 12 spins, as a function of temperature for var- 
ious applied transverse field strengths, where again the 
initial state of the evolution is the canonical equilibrium 
state at a temperature given by /3 = 20. 



Magnetization: The situation for the transverse mag- 
netization remains the same as in the other dimensions: 
While the canonical equilibrium state has zero magneti- 
zation for all temperatures, the time-evolved state, when 
averaged over a large time (away from the initial time), 
produces nonzero magnetization for all nonzero initial 
fields. 

Correlations: The behavior of the T^^ correlation is 
given m the Fig. HI (left). The equilibrium state T^^ 
(curve with black circles) is of course vanishing at /3 = 0, 
which then decreases to —0.963 at low temperatures. The 
long-time average of T^^ in the time-evolved state (with 
the initial state of the evolution being the canonical equi- 
librium state with /3 = 20) has different values for differ- 
ent initial applied fields. They are represented as hor- 
izontal lines in Fig. [6] (left), for a/ J = 0.2, 0.6, 1.2, 
and 2. The qualitative feature is again the same as in 
the other dimensions: Although T^^ is ergodic for low 
applied magnetic fields, it ceases to be so for high mag- 
netic fields, although the nonergodicity obtained is not 
"strong" . The behavior of the T^^ correlation is similar 
to that of T^^ . The T^^ correlation is ergodic for low 
field strengths, while becoming strongly nonergodic for 
high fields. 

Entanglement: We observe that entanglement is er- 
godic for all values of the applied initial magnetic field. 
In Fig. [S] (right), the long-time average of the nearest 
neighbor entanglement is depicted as horizontal lines for 
a/ J — 0.2,0.6, 1.2 and, 2, and all of them intersect the 
corresponding equilibrium curve. 

Note that we find that more classical correlations are 
violating strong non-ergodicity as we increase the dimen- 
sion of the system. 



VIII. CONCLUSION AND DISCUSSION 

We considered the transverse quantum XY model with 
nearest neighbor interactions in low dimensions. We 
found that all the two-body classical correlations along 
with the magnetization in the system, with a suitably 
chosen magnetic field, cannot equilibriate, and hence re- 
main nonergodic, while the bipartite entanglement in the 
system with the same magnetic field can be ergodic. The 
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thesis is true for the infinite chain, which is exactly diago- 
nalizable. The behavior low dimensional finite systems ~ 
ladder and two-dimensional models with periodic bound- 
ary conditions - mimics the one-dimensional infinite one. 

We are therefore led to the conclusion that although 
quantum correlations have traditionally been understood 
to be constituted out of classical correlations and mag- 
netizations, it does not necessarily imply that statistical 
mechanical properties of quantum correlations will be in- 
herited from those of classical correlations and magneti- 
zations. 

Apart from the above implication of fundamental na- 
ture, the result can potentially have implications in the 
decoherence problems in quantum information process- 
ing tasks [22]. 

Lastly, entanglement has been known to be useful for 



information processing tasks, and in these tasks, it clearly 
outperforms all systems which only possess only classical 
correlations [ll-y]- Recent investigations have revealed 
the importance of entanglement-like concepts in dealing 
with statistical mechanical concepts like quantum phase 
transitions 0], which however, can also be detected 
by using classical correlations The results obtained 
in this paper point out that even statistical mechanical 
phenomena, like ergodicity, can appear in quantum en- 
tanglement independently of classical correlations. 
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